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We quantify the important effect of strong final state interactions in the weak K —* 2ti amplitudes, 
using the measured tt-tt phase shifts with J = and / = 0, 2. The final rescattering of the two pions 
provides a strong enhancement of the A7 = 1/2 amplitude, which so far has been neglected in the 
theoretical predictions of e'/e. This correction increases the Standard Model prediction of e'/e to 
values in good agreement with the experimental measurements. 



It is well known that, at centre-of-mass energies 
around the kaon mass, the strong S~wave tt-tt scattering 
generates a large phase shift difference [5^ ~ ^o) {^k) — 
45° ± 6° between the 7 = and I ~ 2 partial waves 
0. In the usual description of K ^ 2n decays, this 
effect is explicitly taken into account, through the fol- 
lowing decomposition of the relevant isospin amplitudes 
with / = 0, 2: 

Ai = A[K->iTnr)i]=Aie<. (1) 

It has also been suggested [||-|| that final state interac- 
tions (FSI) play an important role in the observed en- 
hancement of the 1 = decay amplitude, A0/A2 « 22.2. 
However, their impact on the direct CP-violating param- 
eter e'/e has never been properly estimated. 

At lowest order in Chiral Perturbation Theory 
(ChPT), O(p^), the decay amplitudes do not contain any 
strong phase: 



Ai 



where Cq = 58 
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(2) 



I a/2 527, with 58 and 
527 unknown chiral couplings, corresponding to the two 
lowest-order AS* = 1 operators in the momentum expan- 
sion, transforming as (8l, Ir) and (27^, Ir), respectively, 
under the chiral group [^,1). The phenomenological de- 
termination of these couplings from K — > 27r, taking the 
measured phase shifts into account through (Q), gives 
|Co| ~ 5.1 and I527I ~ 0.29. 

The above procedure is not quite consistent, because 
the strong phases are put by hand. Those phases origi- 
nate in the final rescattering of the two pions and, there- 
fore, are generated by chiral loops which are of higher 
order in the momentum expansion. Since the strong 
phases are quite large, one should expect large higher- 
order unitarity corrections. The existing one-loop anal- 
yses of if ^ 2tt HJ^,0| show in fact that the pion loop 



diagrams provide an important enhancement of the Aq 
amplitude, implying a sizeable reduction (~ 30%) of the 
fitted |58| value. However, the phase-shift Sq predicted 
by the one-loop calculation is still lower than its mea- 
sured value, which indicates that a further enhancement 
should be expected at higher orders. 

Many attempts have been made to compute the am- 
plitudes Ai from first principles |p|,|6|j9|-[l^ . Although 
those calculations have provided encouraging results, we 
are still far from getting accurate predictions. In many 
approaches FSI are not included in the computational 
framework. This is the case of present lattice calcula- 
tions [^3|, which are only able to compute the one pion 
{t^\'Has=i\K) rnatrix elements, or estimates at leading 
order in the 1 /Nc expansion ||6|,p^-p^ (the phases^ are 
zero at leading order). Other approaches ]lO| , p^ , 15 in- 
clude some FSI effects, but in a rather incomplete (and 
untrusty) way; however, the fact that the values found 
for the ^0 amplitude are larger than in the previous two 
methods is a clear indication of the important role of FSI. 

Unitarity and analyticity constraints permit a resuni- 
mation of the effects due to FSI jl^. Let us consider 
the off-shell amplitudes Ai{s), with s = (fc^^ -I- k^^^ 

y the 

invariant mass squared of the final two-pion state. In 
the elastic region Watson final-state theorem relates the 
imaginary part of the amplitudes Ai{s) to the strong 
phases (5/^g(s) of the partial wave amplitudes TJ^q(s) 
which describe the rescattering of the two pions in the 
final state. We can write a 1-subtracted dispersion re- 
lation for the amplitude Ai{s), which has the following 
Omncs solution [0,181: 



Aii.s) = Af^ {s-ml) nji-s), 



(3) 



fli{s) = exp ■ 



(s - ml) 



dz 
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We have imposed the subtraction at s = m^, where the 
amphtude has a zero. The constant A^p is the first 
derivative of Ai{s) at the subtraction point, i.e. the co- 
efficient of the linear term in (s — ) of the Taylor ex- 
pansion around s = m^; up to small higher-order mass 
corrections, it is given by the leading 0{p^) term in the 
chiral expansion. 

The dispersion integral in eq. (^) has an imaginary part 
and a real part, so that the complete amplitude can be 
otherwise written as 



5**0 (-) SJij^s), 



(4) 



where ^i{s) originates from the real part of the disper- 
sion integral, which is directly related to the strong phase 
shift. The factor 5ft/ (s) has never been taken into account 
in the short-distance calculation of the amplitudes Ai{s). 

In this letter we show how the lacking of such term 
3?/ does affect numerically both the ratio A0/A2 and the 
direct CP violation parameter e'/e. The result that we 
obtain is the needed one to move both quantities towards 
the experimental findings. 

The dispersion relation in eq. is strictly valid only 
in the elastic region. Above the first inelastic threshold 
a coupled-channel analysis will be demanded Here 
we limit ourselves to the very simple numerical analy- 
sis in the elastic region, which at s = m|. should give a 
quite good approximation to the exact result, including 
inelastic effects. 

We have used the simple parametrizations of the strong 
phases Sq{s) given in ref. pO| . They reproduce well the 
experimental data up to the KK threshold. Other recent 
parametrizations of the experimental data [|T] performed 
coupled-channels analyses, which we do not need in this 
simple first step calculation. The 1 = 2 amplitude is 
found to be elastic and well reproduced even up to 1.6 
GeV. We then evaluated numerically the dispersion inte- 
grals up to 1 GeV for / = and 1.6 GeV for 1 = 2. 



Taking s 



we find that the on-shell K ^ 2-k 



'if' 



decay amplitudes get the dispersive correction factors: 



5ftn = 1.41 ±0.06 



5ft2 = 0.92 ±0.02. 



(5) 



The central values have been obtained by using the 
parametrizations of (5q(s) quoted as best fits in ref. [ pO[ . 
To estimate the errors we have redone the numerical anal- 
ysis with two other sets of parametrizations, given in 
ref. [ pO[ , which underestimate and overestimate, respec- 
tively, the phase shift data. 

Since we have not included inelastic channels, we have 
not taken into account contributions from energies above 
1 GeV (1.6 GeV for / = 2). This is a quite conservative 
attitude, since those contributions tend to increase (de- 
crease) 5fto {^2) to slightly larger (smaller) values. The 
details of the numerical analysis will be given elsewhere 



The corrections induced by FSI in the moduli of the 
decay amplitudes Ai have resulted in an additional en- 
hancement of the A/ = 1/2 to A/ = 3/2 ratio. 



3fio/5ft2 = 1.53 ±0.07. 



(6) 



This factor multiplies the enhancement already found at 
short distances. This is a quite large correction, not taken 
into account previously, which improves all existing cal- 
culations oi A^p. Taking the 5ft/ correction into account, 
the experimental Ai amplitudes imply the following cor- 
rected values for the lowest-order AS* = 1 chiral cou- 
plings in d): 



ICol «3.6 
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0.32. 



(7) 



This "experimental" numbers are not very far from the 
short-distance estimates obtained in the first of refs. [1). 

It is interesting to compare our results (^ with the 
ones obtained by using the lowest-order estimate of the 
strong phases (5q(s) in ChPT. We performed one step to- 
wards the unitarization of the amplitudes Ai{s) by using 
{U - 93 MeV) 



tan (5q'^(s) 



Ami ( 2s 
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which gives the usual representation of (5o'^(s) at thresh- 
old (s ^ 4m^), where tanf^p'^ ^ Jq'^. For / = the 
lowest-order ChPT parametrization of eq. (||) integrated 
up to 1 GeV gives 5fto = 1.21, sensitively underestimating 
the dispersion integral over the fitted phase as expected. 
For 1 = 2 the situation is similar; lowest-order ChPT 
provides a good description of the strong phase only at 
low energies. We obtain in this case 5ft2 = 0.83, also 
lower than the value in eq. (|^) obtained using the fitted 
phase. Note, however, that the ratio 3fio/5ft2 = 1.46, al- 
though slightly lower than the value (H), stays within the 
quoted error bar. This is in agreement with the observa- 
tion made in ref. |jl[ that the difference 6q — 6q gets smaller 
chiral corrections than the individual phase shifts. 

Implications for e'/e. — 

The most striking consequence of the correction factors 
5fto.2 is a sizeable modification of the numerical short- 
distance estimates for the direct CP-violation parameter 
e'/e. A handful way of writing this quantity, used in all 
theoretical short-distance calculations up to date, can be 
as follows [ pl 



e 



_ p(3/2) 



(9) 



and 



where At = Vj*^^^, the phase $ = $j 
the quantities pf^/^) and pf^/^) contain the contribu- 
tions from the hadronic matrix elements of four-quark 
operators with AI = 1/2 and 3/2 respectively: 



2 



p(l/2) ^ ^ Y.y,{^,) (Q,(^))o (1 - r!,+,0 , 
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Here, (Q,)/ = {{mr)i\Qi\K) , r and a; are given by 



Re^o 



Re^2 
Re^o ' 



and the parameter 



1 (W2)lb. 



(10) 



(11) 



(12) 



parametrizes isospin breaking corrections. It is usually 
set to ^ri+ri' ~ 0.25 with large uncertainties. 

The Wilson coefficient factors ?/i(/i) have been com- 
puted to next-to- leading order accuracy |]lT| , |l3| . Since 
the hadronic matrix elements are quite uncertain theo- 
retically, the CP-conserving amplitudes KeAj, and thus 
the factors r and ui, are set to their experimentally deter- 
mined values; this automatically includes the FSI effect. 
All the rest in the numerator is theoretically predicted 
via short-distance calculations, because the leading con- 
tributions come from the operators Qq and Qg whose 
matrix elements cannot be directly measured from the 
K ^ 2tt decay rates. 

As a consequence, since the relevant matrix elements 
(Qe.s)/ are usually taken from lattice calculations or 
large-iVc estimates , which do not include FSI correc- 
tions, this procedure produces a mismatch with the FSI 
included phenomenologically in the values of r and lo. 
This can be easily corrected, introducing in the numer- 
ator the dispersion factors 'Sti that we have estimated. 
This implies a large enhancement of the predicted value 
of e'/e by roughly a factor of 2. 

A fast way to estimate the numerical enhancement is 
through the approximate formula [pl| 



(1/2) 



ri+n' 



0.4 B, 



(3/2) 



(13) 



where i?g^^^' and Bg^^^'^ parametrize the matrix ele- 
ments of the QCD penguin operator Qg and the elec- 
troweak penguin operator Qg, respectively, in units of 
their vacuum insertion approximation. These parame- 
ters are usually taken to be (from Lattice calculations 
PI and 1/Nc considerations |ll|) S^^^^^ = 1.0 ± 0.3 
and Bg = 0.8±0.2. Since those estimates do not 
include the FSI effect, their values should be multiplied 
by the appropriate factors % and 3?2, respectively [ p3| . 
However, the term i3g^^^''r2,,+,,' in eq. ( p^ ) should be 
multiplied by 5R2 and not by SRq, because it corresponds 
to two final pions with 1 = 2. 

The three terms of eq. (p^ should then be corrected 

to B^'/'Vs/ 



S^^/^^f7^+^/|Fs/ = 0.23 ±0.07. In the latter, we have 
set equal to 0.25, disregarding its (probably large) 

error. This amounts to a FSI enhancement of e'/e by a 
factor 2.1. 

Thus, the so-called "central" value in refs. ||ll|, e'/e = 
7.0 X 10~^ gets increased by the FSI correction to e'/e = 
15 X 10~^, which compares well with the present experi- 
mental world average |Q 



Re (e'/e) = (21.2 ±4.6) x 10" 



(14) 



The final theoretical prediction of e'/e depends on 
other hadronic and quark mixing parameters, which in- 
troduce a rather large uncertainty. It has been said in 
refs. [0 and that the theoretical prediction can only 
be made consistent with the experimental value ( p^ for 
rather extreme values of the input parameters, requiring 
a conspiracy of several inputs in the same direction. The 
important point we want to stress here is that, once the 
effect of FSI has been properly taken into account, the 
experimental value of e'/e can be easily obtained with 
reasonable values |2^ of the different inputs. 

Our calculation of the unitarity 5R/ correction factors, 
provides a hint on why the Trieste and Dortmund 
p5[ groups were able to obtain larger values of e'/e than 
the Munich ||ll| and Rome |jl^ collaborations. The first 
two groups did include some FSI corrections, at the one- 
loop level, which pushed their predictions in the right 
direction. The use of the dispersive factors 3?/ offers a 
more powerful and consistent way to estimate these FSI 
effects, because the Omnes exponential 57/ (s) sums the 
unitarity logarithms to all orders in the chiral expansion. 

More work is still needed in order to get an accurate 
prediction for e'/e. In the meanwhile, our calculation 
demonstrates that, within the present uncertainties, the 
Standard Model theoretical estimates can easily accom- 
modate the experimental value, without any need to in- 
vocate a new physics source of CP violation phenomena. 
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1.4 ±0.3, B^^^^^'Ifsi = 0.7 ±0.2 and 
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